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Abstract. We present Ha line-profile analyses for the seven SU UMa type dwarf novae AK Cnc, WX Get, AQ Eri, VW Hyi, 
RZ Leo, TU Men, and HS Vir All data sets are treated in the same manner, applying a sequence of techniques for each system. 
The basic ingredients of this sequence are the diagnostic diagram to determine the zero point of the orbital phase, and Doppler 
tomography to visualise the emission distribution. We furthermore introduce a new qualitative way of to evaluate the Doppler 
fit, by comparing the line profile of the reconstructed with the original spectrum in the form of the V/R plot. We present the 
results of the analysis in the compact form of an atlas, allowing a direct comparison of the emission distribution in our targets. 
Although most of the data sets were not taken with the intention of a line-profile analysis, we obtain significant results and are 
able to indicate the type of the additional emission in these systems. Our objects should have in principle very similar physical 
properties, i.e. they cover only a small range in orbital periods, mass ratios, and mass-transfer rates. Nevertheless, we find a 
large variety of phenomena both with respect to the individual systems and also within individual data sets of the same object. 
This includes 'canonical' additional emission components from the secondary star and the bright spot, but also emission from 
the leading side of the accretion disc. 

Key words. Accretion, accretion disks - Line: profiles - Methods: data analysis - Astronomical data bases: miscellaneous - 
Novae, cataclysmic variables 



1. Introduction 

Cataclysmic variables (CVs) are close interacting binaries with 
a white dwarf as primary and a late-type main-sequence star 
as secondary component. The latter fills its critical Roche lobe 
and thus enables mass transfer through the inner Lagrangian 
point into the Roche volume of the primary. In the absence of 
strong magnetic fields, the gas stream from the secondary star 
dissipates and forms an accretion disc around the white dwarf 
over which the mass transfer takes place. 

The accretion disc is the dominant light source in these sys- 
tems. It is also the primary origin of the strong emission lines 
especially of the H and He series, which are characteristic for 
a CV spectrum (e.g. Smak ll98l1 Home & Marsh . 1986i Home 
I1995> . In most CVs, the disc is too bright to allow the detec- 
tion of the absorption features of the primary and the secondary. 
The emission lines thus represent the major source of informa- 
tion on the orbital motion. Furthermore, the assumption of a 
Keplerian velocity law allows the derivation of several system 
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parameters like masses and inclination. These parameters are, 
of course, vital for an understanding of the physics in these 
systems, and furthermore represent the only means to test the- 
oretical models on CV evolution (see, e.g. Warner [19951 for a 
comprehensive overview on CVs). 

This approach is only feasible if the emission lines origi- 
nate in a region moving with the same mean velocity vector 
as the primary star. Only in this case do they track its orbital 
motion. However, several examples show that, additionally to 
the disc-borne emission, other emission components distort the 
originally symmetric line profile (e.g. Stover 1981 1. The most 
prominent sources of this emission are the region around the 
bright spot (the place where the gas stream impacts on the disc), 
the gas stream from the secondary star, or an irradiated part on 
the surface of the latter, but also apparently enhanced emis- 
sion regions in the accretion disc, which, in a number of cases, 
make a physical explanation difficult (Tappert 1 999 , North et 
al. l2001> . All these regions have in common that they are 'iso- 
lated' in the sense that they are not centred symmetrically on 
the white dwarf. Consequently, they have both different veloc- 
ity amplitudes and phasing with respect to the primary. The ra- 
dial velocity parameters of the combined line profile as a whole 
thus will not reflect the motion of the white dwarf, leading to 
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erroneous system parameters. We will call these asymmetric 
emission regions in the following 'isolated emission sources' 
(lES). 

Fortunately, in most cases, the distorting components are 
confined to the line centre. ' A careful examination of the line 
profile is thus usually able to separate the influences of addi- 
tional components and allows the determination of the system 
parameters with some confidence. Furthermore, a recent sta- 
tistical study showed that the additional components represent 
more than just a mere obstacle in the derivation of the system 
parameters, but in principle provide valuable information on 
the physical processes in CVs (Tappert & Hanuschik 2001 1. 
Last, not least, the discovery of spiral shocks in CV discs by 
the means of Doppler tomography (Steeghs et al. ll997> proves 
the importance of line profile analyses for the understanding of 
accretion physics. 

Kaitchuck et al. <1994> were the first to re-analyse a large 
sample of 'old', i.e. already elsewhere published, CV spectra 
with a newly developed technique. Their atlas of Doppler maps 
allowed a direct comparison of the various types of emission 
distribution. The atlas included 18 systems, 3 of them with an 
orbital period shorter than 3 hours. 

In a similar spirit, the scope of the present paper is to apply 
a sequence of techniques of line profile analysis to a sample of 
CVs. Our concentration on short-period dwarf novae allows us 
to study the variety of phenomena presented by systems which 
can be assumed to have very similar physical properties, i.e. 
a low mass ratio q - M2/M1 < 0.3 and low mass-transfer 
rates of the order M ~ lO^'^Mo/yr. Furthermore, according 
to standard C V evolution theory, low mass transfer rates imply 
optically thin discs. These are strong line emitters and could 
reveal structural details more easily than optically thick discs. 

The present atlas also provides a significant increase of 
Doppler maps for short-period CVs, whose number of avail- 
able line-profile analyses is in strong mismatch with respect to 
the number of known systems, compared with long-period CVs 
(Tappert & Hanuschik 200 U . This is basically due to the fact 
that for these faint systems the orbital period is easily accessi- 
ble for small telescopes through the superhump phenomenon, 
while detailed studies need large telescopes. The latter there- 
fore usually concentrate on a few 'interesting' objects, like WZ 
Sge (e.g., Skidmore et al. 2000 and references therein). Our 
present sample instead is not part of this 'selection effect', as 
none of the systems has (yet) been the subject of a more de- 
tailed investigation involving the necessary large telescopes. It 
therefore adds to a more 'unbiased' view of the properties of 
SU UMa stars. 

Following our goal of direct comparison, we restrict our 
analysis to the Ha profile, which also represents the most ac- 
cessible emission line. 

2. The sample 

The systems in question are the seven SU UMa type dwarf no- 
vae AK Cnc, WX Cet, AQ Eri, VW Hyi, RZ Leo, TU Men, 

' The high velocity wings in SW Sex systems are an important ex- 
ception (e.g. Hoard et al. .2000» and references therein). 



and HS Vir With the exception of VW Hyi, the data were not 
taken with the intention to perform a line-profile analysis, but 
to determine the orbital period. The observational background 
and the properties of the data have been published elsewhere, 
together with the results of the initial analysis. In Tabled we 
therefore present only a summary of the basic parameters of 
the data sets, and refer the reader to the respective publications 
listed in the table for further details. 

The available data on HS Vir actually consists of more than 
the one set covered here. However, the system showed large, 
apparently erratic, variations from one set to another, and also 
within individual data sets. A more detailed analysis is required 
in this case, which is beyond the scope of the present paper, 
since it aims at a uniform treatment of the data. We have there- 
fore included only the 'most well behaved' data set of this ob- 
ject as an example. 

3. Analysis sequence 

The following sequence of analysis techniques is applied to 
each data set: 

1. Visual inspection: Apart from the changes due to the orbital 
motion and/or to additional emission sources, the emission-line 
profile is also aff'ected by luminosity changes of the system. 
Roughly, the strength of the disc-borne emission lines is ex- 
pected to decrease with increasing disc brightness, due to the 
increasing optical thickness of the accretion disc. Additional 
emission components, however, might be affected quite dif- 
ferently by such a variation, so that their influence becomes 
stronger (e.g., emission from the secondary star) or weaker 
(e.g., emission from the bright spot). Line profiles of the same 
CV in different states thus usually show quite different phe- 
nomena (compare, e.g., Martfnez-Pais et al. ll994lll996l for SS 
Cyg in quiescence and outburst, respectively). A visual com- 
parison of line profiles of different data sets therefore provides 
a first hint as to whether a combination of these data is possible 
in order to increase the phase resolution, or if each set has to be 
treated individually. 

Similarly, in high inclination systems, the brightness 
variation due to the orbital motion can affect the line profile 
at specific orbital phases. This is especially true for eclipsing 
systems, where also the line emission source is partly or 
completely obscured (e.g. Kaitchuck et al. '1998't. These parts 
of the data sets usually have to be excluded from the line 
profile analysis. 

2. Diagnostic diagram: This method makes use of the double- 
Gaussian convolution with the line profile (Schneider & Young 
[T980 Shafter lT983bl . which allows to measure radial velocities 
of different parts of the line. These velocities v, (1^) can be fitted 
with a sinusoidal function 

Vr(<p) = r - ■'^em sin(27r(vP - (fio)), (1) 

to yield the systemic velocity y, the semi-amplitude /Tem, and 
the fiducial phase corresponding to superior conjunction of the 
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Table 1. General properties of the data. The dates of observations and the corresponding data sets (numbers usually refer to 
individual nights; for exceptions see text) are presented in Col. 2 and 3, respectively. Column 4 gives the total number of spectra 
available for the analysis. Col. 5 their spectral resolution (FWHM), Col. 6 the telescope used for the observations, and the last 
column lists references concerning the previous publication of the data sets. 



star 


dates 


sets 




AA [A] 


telescope 


reference 


AK Cnc 


1995-03-19, 21, 22 


1-3 


40 


4.5 


2.5m LCO 


Arenas & Mennickent ( 1998 1 


WX Get 


1992-12-07, 08, 10 


1-3 


37 


2.0 


2.2m ESO/MPI 


Mennickent 1 1994j 


AQEri 


1992-12-11 


l(a-c) 


44 


2.0 


2.2m ESO/MPI 


Mennickent ( 1995^} 


VW Hyi 


1994-10-15 to 17 


1-3 


64 


1.1 


1.9m SAAO 


Tappert 1 19991 


RZLeo 


1995-02-10, 11, 03-21, 23 


la, lb, 2a, 2b 


46 


6.0, 4.5 


2.2m ESO/MPI, 2.5m LCO 


Mennickent & Tappert <2001> 


TUMen 


1992-12-09, 10, 12, 14 


1-4 


67 


2.0 


2.2m ESO/MPI 


Mennickent 1 1995a} 


HS Vir 


1998-05-30 


1 


34 


2.5 


3.5m NTT 


Mennickent et al. <1999> 



primary star (po. When these parameters are plotted as a func- 
tion of the separation d of the two Gaussians, a so-called di- 
agnostic diagram is obtained, which illustrates the variation 
of the parameters as different parts of the line are measured 
(Shafter 1983a I. In order to evaluate the significance of the re- 
sulting variation, the corresponding error of the semi-amplitude 
o"(A'em)/^em IS included in the plot. A sharp increase of this 
value at large separations indicates the point beyond which the 
noise begins to dominate the signal from the line wings and 
the parameters become unreliable. The values belonging to the 
largest separation d^ax before this happens are usually chosen 
as the best possible approximation to the correct ones, so that 
^em - Ki, i.e. the true semi-amplitude of the orbital motion of 
the primary. 

A second criterion is based on the expectation that the pa- 
rameters should (again) approach constant values at large sepa- 
rations when additional emission sources cease to influence the 
line profile. This is usually the stronger criterion, as the first 
one is often ambiguous. 

The diagnostic diagram furthermore allows us to fix the 
fiducial phase, which is a prerequisite for the evaluation of the 
subsequent Doppler mapping. 

3. Doppler tomography: 



This technique, introduced by Marsh & Horne (I1988> to 
the realm of interacting binaries, interprets the one-dimensional 
snapshots represented by a spectrum taken at a specific orbital 
phase as projections of a two-dimensional emission distribu- 
tion. The latter is reconstructed as a map in velocity space via 
an inversion of these projections (see Marsh .2001 . for a recent 
review). 

In this work, the implementation of Spruit (1998) has been 
used. We have replaced the original IDL routines by a corre- 
sponding MIDAS interface, but still use the FORTRAN core 
program (version 2.3.1), to run the computation on a Linux 
PC. The input spectra have first been continuum subtracted, 
the extracted line profiles have then been normalised in flux 
(in order to avoid too small numbers for the Doppler routine, 
the total line flux has been set to 100) to minimise artefacts 
due to intensity variations. A smearing kernel of 2 pixels was 
used to compute the Doppler maps, which were subsequently 
smoothed with a 3 x 3 pixel (11 x 11 in the case of V W Hyi) 



averaging filter for noise reduction. Tests with larger kernels 
showed that the visual information content of the resulting map 
basically remains the same, but that the profile comparison (see 
step 4) suffered from too strong smoothing. A kernel size of 1 
pixel, on the other hand, in almost all cases led to infinite itera- 
tion loops due to the noise becoming dominant. All other input 
parameters for the Doppler routine were set to default values 
(SDruit .l9ggl . 

The y velocity resulting from the diagnostic diagram was 
taken as a first guess for correcting the input line to its rest 
wavelength. It was then iteratively adjusted by comparing the 
reconstructed and original average spectra. 

The resulting Doppler maps are orientated as usual, i.e. 
the fiducial phase is at = and v,, > and increases in 
clockwise direction. If the phase corresponds to the orbital 
motion this results in the white dwarf being located at (v^,v^.) 
- {0,-Ki) and the secondary star at (0,7^2). 

4. V/R-plot: Perhaps the simplest way to examine a line pro- 
file for asymmetry is to compare its blue part to its red one. 
This has been done in the past for double-peaked profiles by 
computing the ratio between the intensity of the blue (violet) 
peak and the one of the red peak, the so-called V/R ratio (e.g., 
Mennickent I 1 9941 1 . In order to apply this method also to single- 
peaked profiles, one can define a different, but equivalent, V/R 
as (Tappert [D^ 



V/R = log 



F{V) 
F(R) ' 



(2) 



where F(V) is not the peak intensity but the flux of the blue 
part and F(R) the corresponding one of the red part. The point 
which separates both halves, A^, is chosen as the centre of the 
line flanks at a specific intensity value which has to be low 
enough to avoid being affected by the central valley in double- 
peaked profiles, and high enough to avoid the noise-dominated 
extreme line wings. The sigma of the individual V/R value is 
computed by a Monte Carlo simulation which adds a random 
value, uniformly distributed within an interval determined by 
the S/N, to each data point in the spectrum, and measures the 
V/R for a thousand of such newly computed data sets. 

We here use the V/R plot to compare the original spectrum 
with the reconstructed one from the Doppler map. In this way it 
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Table 2. Radial-velocity parameters for Ha as derived by the 
analysis of the diagnostic diagram. Col. 2 gives the FWHM of 
the Gaussians used for the diagram. Col. 3 the critical sepa- 
ration. Col. 4 and 5 contain the resulting semi-amplitude and 
the error of the zero phase, respectively. The y parameter is 
not listed, as it reflects more the wavelength calibration pro- 
cess than a real physical parameter. 



data set 


FWHM 






cr(<po) 


Fig. 




[A] 


[A] 


[km/s] 


[orbits] 




AK Cnc 


4 


34 


42(09) 


0.036 




WXCet 


4 


52 


59(18) 


0.061 


1X31 


AQEri 


2 


44 


36(09) 


0.042 


lAJl 


VWHyi 


4 


40 


38(09) 


0.039 


lAin 


RZLeo 1 


4 


66 


49(13) 


0.042 




RZLeo 2 


4 


62 


51(15) 


0.050 


lAlSl 


TUMen 


2 


48 


84(09) 


0.017 


lAJSl 


HS Vir 


2 


64 


98(09) 


0.017 





can be quantitatively evaluated if the Doppler tomography has 
been able to reproduce the general shape of the individual line 
profiles. 

In the atlas presented in Appendix IaI we furthermore include 
two-dimensional greyscale plots of the original and the dif- 
ference data, the latter being computed by subtracting the re- 
constructed spectrum from the original one. In this way, also 
possible intensity deficits (black features in the plot) or ex- 
cesses (white features) of the reconstructed data become vi- 
sually clear. 



4. Results 

4.1. Splitting of data sets 

For all objects, the analysis methods were first conducted on 
each individual data set (i.e., corresponding to one night of ob- 
servation; see Table 0. Where no significant diff'erences were 
found, the analysis was repeated on the combined set in order 
to minimise noise. 

The data of AQ Eri represent a special case. Here, the line 
profiles proved to vary both in strength and in shape on non- 
orbital timescales (compare especially phases 0.1 and 0.8 in 
the upper plot of Fig. lA.7> . We therefore split the set in three 
parts la, lb, Ic, with each set covering one full orbit. Set lb 
overlaps with the other two (it includes 4 spectra from set la 
and 3 from set Ic), but sets la and Ic are independent of each 
other 

The appearance of the line profiles of HS Vir (upper plot 
of Fig. IA.20> suggests a similar treatment. However, here the 
phase distribution of the data did not allow a division into in- 
dependent sets which would cover a complete orbit. We con- 
ducted the analyses for different subsets, but did not find any 
significant differences with respect to results or residuals. 



4.2. Diagnostic diagrams 

The diagnostic diagrams proved to be robust with respect to the 
differences between the individual data sets. With the exception 
of RZ Leo, where the insufficient precision of the orbital period 
does not allow for an interpolation of the phases, we therefore 
used the combined data sets. 

In the present application of the diagnostic diagram we are 
mostly interested in determining the zero point of the orbital 
phase, on which depends the interpretation of the subsequent 
Doppler maps. For the majority of the objects in our sample, 
this choice was unambiguous, as after clearing distortions due 
to structure in the line core this parameter remained constant 
over a large range of Gaussian separations. We therefore based 
our choice of the maximum separation t/^ax on the stability of 
the other parameters (of y in the case of AK Cnc and WX Cet, 
y and K^^ for TU Men) where the noise parameter a-{Ksia)/Kem 
appeared not to be reliable, on a combination of the noise and 
the constancy criteria (AQ Eri), or (for HS Vir) only on the 
noise behaviour. In the cases of VW Hyi and RZ Leo, ipo does 
not reach constancy within the allowed range with respect to 
the noise criterion. We here chose the values corresponding to 
the largest separation within this range. Judging from the ap- 
pearance of the variation of (po at the separations before this 
point (slight monotonic decline), we do not expect the error in 
the determination of the fiducial phase to be larger than ~0. 1 
orbits. 

For completeness, we have listed the parameters corre- 
sponding to our choice of t/^ax in Tabled 

4.3. Doppler maps 

The general appearance of the Doppler maps is a rather noisy 
one. This is due to the fact that the S/N of the data could not be 
improved in the majority of our systems by co-adding phases 
of different orbits because of the observed short- and long-term 
variations. For the same reason also the phase coverage often 
stays incomplete and/or non-uniform. All this gives rise to arte- 
facts in the Doppler maps, and one has therefore to be careful 
not to over-interpret certain features. 

Nevertheless, the tests by the means of the trailed spec- 
trograms and the V/R plots show that the reconstructed spec- 
tra produced by the Doppler fitting agree fairly well with the 
original data, with the exception of HS Vir (see Sect. I4.3.7> . 
Differences in the V/R values (e.g., lower plots for WX Cet, 
Figs. lA.4l and lA.6l RZ Leo, Fig. lA. 171 can largely be attributed 
to noise in the original spectra and to the fact that due to the 
Doppler smoothing (or smearing) the S/N in the reconstructed 
spectra is often up to 5 times higher. We therefore conclude 
that the line profiles themselves have been well reproduced. As 
shown by the trailed spectrograms, in general the same can be 
said for the intensity behaviour, with a few notable exceptions 
which are discussed below for the respective systems. 

4.3.1. AKCnc 

The Doppler map shows a symmetric emission distribution 
(Fig. lA.2> . AK Cnc therefore appears to be one of the few sys- 
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Fig. 1. Example of an original continuum-normalised spectrum 
of the WX Cet data set (bottom) and its corresponding model 
spectrum (top). The latter has been displaced vertically by 2.5 
intensity units. 



tems which show no presence of any lES; the survey of Tappert 
& Hanuschik (2001 ) found only 4 such objects in a sample of 
68 CVs (not including AK Cnc). Unfortunately, the system was 
not observed in quiescence, but in decline from outburst. It thus 
remains unclear if the absence of IBS is caused by a still opti- 
cally thick disc (i.e. if this is only a temporary absence), or 
simply due to the low inclination of the system. 

4.3.2. WX Get 

The Doppler maps for the WX Cet are amongst those with the 
most noisy appearance (Figs. lA.4IA.6l l. The low S/N of the data 
is probably also responsible for the fact that the clear signature 
of disc emission, i.e. a ring-shaped emission distribution cen- 
tred on the primary, is not seen, although the line profiles are 
double-peaked (Fig. lA.3t . Instead all maps show a rather spotty 
emission distribution, however, with a maximum on the leading 
side (i.e. Vj > 0) as a common feature. 

The low S/N and the non-uniform phase distribution make 
the WX Cet data sets the ideal test case to examine the signif- 
icance of the features in the Doppler map. For this, we estab- 
lished a set of model spectra using the line profile prescription 
by Smak ( II981> . We adjusted the model parameters until a vi- 
sually good agreement with the shape and quality of the origi- 
nal WX Cet data was reached. The model data were displaced 
in wavelength according to Eq. Q, using the phases of the WX 
Cet data and the parameters obtained by the application of the 
diagnostic diagram {Ki - 59 km s"') and the Doppler tomogra- 
phy (y - 260 km s"'). Fig. Ogives an example of the resulting 
model spectrum. 

These spectra represent a symmetric disc emission, without 
any additional component. The corresponding Doppler maps 
(Fig.O, however, show the same spotty appearance as for the 
original data set, with the ring-shaped disc signature being only 
slightly more evident. The low intensity arcs and blobs in many 
of the Doppler maps (see also the VW Hyi data, especially Fig. 
IA.13> can also be explained by the data quahty. It even ap- 
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Fig. 2. Doppler maps for the model spectra with respect to the 
WX Cet data sets 1-3 (top to bottom). 



pears not improbable that the individual maps of WX Cet could 
be simply reproduced by noisy spectra without the need to in- 
voke an additional component (compare especially the map for 
model set 2 with the WX Cet data). However, the fact that all 
WX Cet maps show an emission maximum on the leading side 
of the disc indicates that this represents a real feature. 
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Fig. 3. The Ha continuum magnitude of the VW Hyi data as 
a function of phase. The upper plot has been composed out of 
data sets 1 (o), 2 (a), and 3 (□), with the respective average 
values subtracted. The lower plot gives those three data sets 
averaged into bins of 0.05 orbital phases, with the error bars 
referring to +lcr deviation of the averaging process. The data 
are repeated to show two orbits. 



4.3.3. AQ Eri 

The maps for the individual data sets show significant dififer- 
ences (Figs. IA.8IA.10l . with an emission maximum close to 
the possible position of the secondary star at (0, K2) moving 
to a position in the lower left quadrant of the Doppler map. 
However, the important difference to the behaviour in other 
systems (e.g. VW Hyi) is that this evolution takes place over 
~2.5 orbits, i.e. on very short timescales. It is therefore doubt- 
ful if the maps in this case are representative of the emission 
distribution at all, or if the timescales of the profile variations 
in AQ Eri are shorter than even one orbit. The fact that the com- 
parison between the original and the reconstructed data shows 
evidence of systematic differences (especially the grey scale 
plot for set lb and the V/R plot for set Ic), seems to indicate 
the latter. 



4.3.4. VW Hyi 

The maps show an evolution of the emission distribution from a 
combined secondary star / gas stream location in set 1 to a dom- 
inating bright spot in set 3 (Figs. lX"l2IA.14> . The ring-shaped 
disc signature is present, but - like in WX Cet - distorted due 
to the low S/N. 

While the V/R plot comparison of the original and the re- 
constructed data are in good agreement with each other, there 
are some systematic differences to be seen in the trailed spec- 
trograms especially for set 1. This concerns a narrow emis- 
sion component that is stronger in the reconstructed data from 
phases 0.7 to 1.2 and weaker during the other phases. We at- 
tribute this to an occultation effect of this emission component, 
the latter coming most probably from the gas stream or the sec- 
ondary star (which shows us its non-irradiated backside around 
phase 0.0). This violates one of the prerequisites for Doppler 
mapping, specifically that all emission sources have to be vis- 
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Fig. 4. Long-term light curve for VW Hyi based on AAVSO 
observations (Mattel 119991 private communication). Upper 
brightness limits are marked by arrows. The dates of the spec- 
troscopic measurements are indicated by the dashed lines. 



ible at all phases. We find a similar phenomenon for TU Men 
(Section R3^ . 

While we did not detect any systematic variation in the Ha 
equivalent widths, the continuum - i.e. the optically thick emis- 
sion-shows a hump at phases 0.75- 1.0 -probably from the gas 
stream -, followed by a shallow eclipse up to phase 1.2 (Fig. 
^} Although the behaviour of the continuum may be differ- 
ent from that of the emission line, this indicates that the in- 
clination of the system is high enough for occultation effects. 
Also the fact that the differences in the trailed spectrograms de- 
crease from set 1 to set 3, while at the same time the Doppler 
maps show that emission from the secondary star becomes less 
dominant in this sequence, supports the interpretation that the 
narrow emission component originates from the secondary star 
and is occulted during the corresponding phases. 

The data presented here have been taken 21 days after max- 
imum of a normal outburst, with the system apparently being 
in a long, slow decline (Fig.|4}. The phenomena reported here 
might therefore either represent a continuous evolution from 
the outburst to the quiescent state, or a short-term fluctuation 
(timescale of a few days). 

4.3.5. RZ Leo 

Both data sets yield very similar emission distributions in the 
form of two maxima on the leading and on the trailing side of 
the disc (Figs. IaT61 and IaTTtI . with the latter maximum being 
the stronger one. The intensity difference appears to be more 
pronounced in data set 2, but it remains unclear if this variation 
is significant or due to noise. 

In view of the spotty appearance of other Doppler maps in 
this atlas (WX Cet, VW Hyi) it is justified to ask if the spots 
in RZ Leo might be artefacts as well. However, in this case. 



^ The uncertainties in the determination of the continuum (which 
was measured with a graphics cursor) are hard to quantify, but the fact 
that all three data sets show the same basic features at identical phases 
proves the reality of this variation. 
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also the trailed spectrograms show the unambiguous presence 
of an additional component. At least the spectrogram for data 
set 1 furthermore shows evidence in the form of a double wave 
for a second component. While the second spot (on the leading 
side) in data set 2 might therefore be questionable, in set 1 the 
strength of this feature appears to exclude an artificial origin. 

4.3.6. TU Men 

The Doppler map of TU Men shows a strong emission feature 
from the secondary star and otherwise pure disc emission (Fig. 
IA.19> . The differences between the reconstructed and the orig- 
inal data which can be seen in the grey scale plots affect the 
same phases as in the case of VW Hyi (Section |4.3.4> . Also 
the explanation is the same: if an emission feature is not vis- 
ible at all phases (i.e. in violation of one of the prerequisites 
for Doppler tomography), the Doppler fit will smear the signal 
over all phases. The feature will therefore be artificially en- 
hanced during the phases when it is absent in the original data 
(here: phases 0.8-1.2), and it will be reconstructed with dimin- 
ished strength in the phases where it is visible (here: 0.3-0.7).^ 
However, while for TU Men the residuals are more posi- 
tive than negative, the VW Hyi data set 1 shows the opposite 
behaviour In other words, in TU Men the IBS feature is dimin- 
ished rather than enhanced, while in VW Hyi the enhancement 
is stronger. 

A possible reason could be a much stronger obscuration in 
VW Hyi. Indeed the differences in the strength of the emission 
feature at the respective phases appear to be more pronounced 
in VW Hyi than in TU Men. Unfortunately, the TU Men data 
had not been flux-calibrated so that the strength of a possible 
obscuration effect in TU Men cannot be compared to that in 
VW Hyi. 

4.3.7. HS Vir 

The Doppler fitting resulted in a dominant emission feature in 
the upper left quadrant of the map (Fig. lA. 211 . However, both in 
strength and in shape there are significant residuals when com- 
paring the reconstructed with the original data. Possible rea- 
sons are either an obscuration effect like in VW Hyi (Section 
I4.3.4> or sub-orbital variations like probably in AQ Eri (Section 
I4.3.3> . A comparison of the line profiles of the individual phase 
bins (Fig. IA.20t makes the second possibility appear more 
probable. A good example is phase 0.15, where the line pro- 
file over two orbits (left and right plot) has changed both in its 
general shape and even in the width of its wings. Furthermore, 
HoweU et al. <1990> found strong variations on non-orbital time 
scales in the light curve of HS Vir, that could not be explained 
by flickering. In this case, Doppler mapping with respect to 
the orbital period yields no representative image of the emis- 
sion distribution. It might be possible that this technique can 
be applied successfully, if a very large data set is used, so that 
non-orbital variations are averaged out. 

^ An extension of the Doppler tomography method to account for 
emission components with variable strength has been very recently 
proposed by Steeghs t20(i?l . 



5. Discussion 

The systems analysed in this paper have in common that they 
are all of the SU UMa subclass, i.e. they are dwarf novae with 
orbital periods (in this case) ranging between 1 .4 and 2.8 h, and 
they have low mass-transfer rates resulting in comparatively 
faint, optically thin, accretion discs. Excluding the extremely 
long periodic system TU Men, the period range spanned by the 
systems even amounts to only 27 min. It is therefore rather sur- 
prising that in these 7 systems with supposedly common basic 
properties we are confronted with almost the complete variety 
of isolated emission sources (lES). 

One can now attempt to look for specific differences be- 
tween the systems regarding their parameters. In general, there 
appears to be a slight dependence for the lES on the inclina- 
tion / of the CV, in that emission from the trailing side of the 
disc (v'l < 0; bright spot, etc.) is more likely to be present in 
high-inclination systems, and emission from the leading side is 
more frequent in low-inclination objects (Tappert & Hanuschik 
1200 U . Other possibly interesting parameters are the mass ratio 
q - M2/M1, which contains information on the stellar compo- 
nents, and the outburst recurrence time free, which is related to 
the transfer and accretion process. Most of these parameters di- 
rectly or indirectly affect the equivalent width Whu of the line. 
If the same mixture of effects is also responsible for the emis- 
sion distribution, one should see a corresponding correlation. 
Finally, the orbital period Poib is indicative of the mass of the 
secondary star 

We have listed these parameters together with the detected 
lES types in Table |3l With the range in f orb being rather nar- 
row, we do not expect it to play an important role in the lES 
distribution for the cuiTent sample. Instead the equivalent width 
seems to qualify as a better sorting criterion. It is influenced by 
a mixture of system (e.g. inclination) and disc (e.g. accretion 
rate) dependent parameters, which is Ukely to be also the case 
for the emission distribution. As all data sets for a specific sys- 
tem show in general similar additional emission components, 
we have sorted the data in Table |3] according to the average 
equivalent width of a system. 

Still, there does not appear to exist any clear correlation 
with the lES type. For example, emission from the trailing side 
(v I < 0) is observed almost over the whole range of Wua, while 
emission from the secondary star (TU Men) is right in the mid- 
dle of our sample. However, we note that emission from the 
leading side is only present in two systems with Who > 100 
A (RZ Leo & WX Cet). A look at the only other parameter 
that is roughly known for all the systems, the recurrence time, 
seemingly strengthens this tendency, since these two WZ Sge 
related systems represent the only ones that exclusively show 
superoutbursts. However this very probably is an artifical ef- 
fect due to the smallness of our sample. lES on the leading side 
have also been at least suspected** e.g. in SU UMa (Thorstensen 
et al. 19861 and CY UMa (Marti'nez-Pais & Casares |l995}. 



Doppler tomography has been published for none of these sys- 
tems, so that the evidence on their emission distribution is only sec- 
ondary. This emphasises the need for more detailed line-profile studies 
for short period CVs. 
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Table 3. Isolated emission features detected in the individual data sets, which here appear sorted with increasing equivalent 
width. Cols. 2 to 4 give basic parameters of the systems, i.e. the period, the inclination, the mass ratio q - M2/M1, and the 
(normal) outburst recurrence time, except for RZ Leo and WX Cet, which show only superoutbursts. These parameters were 
taken from the Ritter & Kolb (,1998 1 catalogue except where indicated. Cols. 5 and 6 contain the data set specific information 
on the equivalent width of the Ha emission line and the accretion state during the observations (d = decline, q - quiescence), 
respectively. The x-velocities in Cols. 7 and 8 refer to the emission maxima in the respective halves of the Doppler maps, other 
emission regions are denoted by their identification. The presence of a feature is indicated by '+', its absence by uncertain 
ones by '?', and very uncertain ones by '??'. 
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which both show also normal outbursts and have recurrence 
times much less than one year (Ritter & Kolb'1998). 

For the other parameters the list is largely incomplete. This 
is basically due to the fact that the calculation of / and q needs 
some knowledge about the secondary star, the latter being dif- 
ficult to detect in short-period CVs. Furthermore, those param- 
eters are usually rather uncertain with 1 cr errors in the order of 
10 to 20%. 

An additional problem is that we are always dealing with 
a combination of parameters, and a sample of 7 systems is 
certainly too small to provide significant statistics in order to 
separate the individual influences. So, while the variety of the 
phenomena detected here might come as a surprise, the lack of 
a clear correlation with the system parameters probably should 
not. 

Perhaps the only system whose emission distribution can be 
fixed to a parameter is TU Men, which shows a strong contribu- 
tion from the secondary star in its quiescent state. In SU UMa 
systems this phenomenon has been hitherto observed only in 
outburst. It appears that in this respect TU Men already quali- 
fies as a long-periodic system, whose disc is bright enough to 
induce emission on the surface of the secondary star. 

Finally, we would like to point out that at least in one sys- 
tem the variations between different data sets are not irregu- 
lar. The analysis of spectra from subsequent nights of VW Hyi 
show a clear long-term evolution from almost pure gas-stream 
emission to a dominant bright spot. At the time of the observa- 
tions, the system apparently is in a state of a prolonged, slow 
decline (with a rate of roughly 1 mag in 40 days) about 20 days 
after a normal outburst and -35 days before a superoutburst 



(Fig. 1^. The long-term variations in the line-profile therfore 
take place during a period where the overall brightness of the 
system is barely changing. Apparently even variations on small 
scales - in the sense that they are not accompanied or immedi- 
ately followed by large photometric variations (i.e., outbursts) 
- of the mass-transfer rate and/or in the disc itself yield signif- 
icant changes in the line profile. 

6. Conclusions 

The here presented atlas of line-profile analyses for 7 SU UMa- 
type dwarf novae increases the number of available Doppler 
tomographies for such systems by almost 100 %. Our maps 
certainly suff'er from the low quality of the data with respect to 
spectral resolution and low S/N. The reason for the former is 
that the data were not taken with the aim to perform this kind 
of analysis, but rather to derive the most basic parameter for 
CVs, the orbital period. Nevertheless, they still served to gain a 
basic idea on the occurrence of isolated emission sources (lES) 
in these systems. 

The noisy appearance of the maps, on the other hand, is 
mostly due to the fact that several systems showed non-orbital 
variations of the emission distribution, with the consequence 
that the respective data sets could not be combined in order to 
improve the S/N, thus actually representing one of the results 
of our study. Apart from that we can draw the following con- 
clusions: 

1 . There is a large variety of lES phenomena present even in 
systems with in principle very similar physical properties. 
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2. With HS Vir and probably also AQ Eri we find two sys- 
tems which show sub-orbital variations in their line profile. 
Consequently, Doppler tomography does not yield a repre- 
sentative image of the emission distribution in these cases. 
Very large data sets, where this type of variation can be 
averaged out, might be needed for this purpose. These sys- 
tems also represent interesting targets for the recently intro- 
duced approach to map the emission line flickering (Diaz 

3. On the other hand, the non-orbital variations in VW Hyi 
occur on long-term timescales and therefore represent a 
trend in the evolution of the emission distribution. Here, 
long-term monitoring of the inter-outburst state via Doppler 
mapping has the potential to provide important clues on the 
outburst mechanisms (e.g. the role of irradiation). 

4. TU Men can now claim the record of being the C V with the 
shortest orbital period (and the only SU UMa dwarf nova) 
to show emission from the secondary star in quiescence.^ 
Previously, this was the nova-like DW UMa (foib - 3.28 h), 
which showed this type of emission in its low state (Dhillon 
et al. [T994l l. 

5. The re-examination of 'old' data with new techniques is im- 
portant and meriting, even if only to motivate more detailed 
studies on the respective systems. 
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Fig. A.l. Top: Line profiles of selected 10% phase bins for AK 
Cnc. Data sets 1 - 3 are plotted from left to right. The time se- 
quence within one data set is symbolised by the sequence solid 
- dashed. Bottom: Diagnostic diagram for AK Cnc, already 
corrected for the derived zero phase. The dotted lines mark the 
chosen separation and the corresponding parameters. 
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Fig. A.l. Doppler Tomography of AK Cnc. Top: Original (left) 
and difference (original-reconstructed; right) spectrum. The 
intensity bar on the right refers to the latter plot. Middle: 
Doppler map. Contour levels are at 90, 75, and 60 per cent of 
the maximum intensity. Bottom: V/R plot of the original (o) 
and the reconstructed (x and dashed line) data, and the residu- 
als (bottom of the plot). In phases to 1 , error bars are given for 
the original data, phases 1 to 2 show those for the reconstructed 
data. 



Tappert et al.: An atlas of line profile studies for SU UMa stars 



11 



Mm 








M I I I I I M M I I I I I I I I I 








-ana -a o a — a o h 
Facial Vfloelty [10^ km/s] 



S 100 




150 



100 



r-. 0,3 




I 1 1 1 1 I 1 1 1 1 1 1 I 1 1 1 I I 1 1 1 I I 1 1 1 I 1 1 I 



o - 



I 1 1 1 1 I 1 1 1 1 1 1 I 1 1 1 I I 1 1 1 I I 1 1 1 I 1 1 1 1 





0:3 
^ti- 



ll I I I I I I I I I I I I I I I I I I I I I I I I I I I 
10 23 30 40 50 60 70 
Separation 

Fig. A.3. Top: Line profiles of selected 10% phase bins for WX 
Cet. Data sets 1 - 3 are plotted from left to right. The time se- 
quence within one data set is symbolised by the sequence solid 
- dashed. Bottom: Diagnostic diagram for the combined data 
set of WX Cet, already corrected for the derived zero phase. 
The dotted lines mark the chosen separation and the corre- 
sponding parameters. 
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Fig. A.4. Doppler Tomography for WX Cet, data set 1 . Top: 
Original (left) and difference (original-reconstructed; right) 
spectrum. The intensity bar on the right refers to the latter plot. 
Middle: Doppler map. Contour levels are at 86, 77, and 68 per 
cent of the maximum intensity. Bottom: V/R plot of the origi- 
nal (o) and the reconstructed (x and dashed Une) data, and the 
residuals (bottom of the plot). In phases to 1, error bars are 
given for the original data, phases 1 to 2 show those for the 
reconstructed data. 




Fig. A.5. Doppler Tomography for WX Cet, data set 2. Top: 
Original (left) and difference (original-reconstructed; right) 
spectrum. The intensity bar on the right refers to the latter plot. 
Middle: Doppler map. Contour levels are at 90, 81, and 72 per 
cent of the maximum intensity. Bottom: V/R plot of the origi- 
nal (o) and the reconstructed (x and dashed line) data, and the 
residuals (bottom of the plot). In phases to 1, error bars are 
given for the original data, phases 1 to 2 show those for the 
reconstructed data. 



Fig.A.6. Doppler Tomography for WX Cet, data set 3. Top: 
Original (left) and difference (original-reconstructed; right) 
spectrum. The intensity bar on the right refers to the latter plot. 
Middle: Doppler map. Contour levels are at 93, 84, and 75 per 
cent of the maximum intensity. Bottom: VIR plot of the origi- 
nal (o) and the reconstructed (x and dashed Une) data, and the 
residuals (bottom of the plot). In phases to 1, error bars are 
given for the original data, phases 1 to 2 show those for the 
reconstructed data. 
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Fig. A.7. Top: Line profiles of selected 10% phase bins for AQ 
Eri. Data sets la and Ic are plotted on the left and on the right, 
respectively. The time sequence within one data set is symbol- 
ised by the sequence solid - dashed. Bottom: Diagnostic di- 
agram for the combined data set of AQ Eri, already corrected 
for the derived zero phase. The dotted lines mark the chosen 
separation and the corresponding parameters. 



Fig.A.8. Doppler Tomography for AQ Eri, data set la. Top: 
Original (left) and difference (original-reconstructed; right) 
spectrum. The intensity bar on the right refers to the latter plot. 
Middle: Doppler map. Contour levels are at 90, 84, and 78 per 
cent of the maximum intensity. Bottom: V/R plot of the origi- 
nal (o) and the reconstructed (x and dashed Une) data, and the 
residuals (bottom of the plot). In phases to 1, error bars are 
given for the original data, phases 1 to 2 show those for the 
reconstructed data. 
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Fig.A.9. Doppler Tomography for AQ Eri, data set lb. Top: 

Original (left) and difference (original-reconstructed; right) 
spectrum. The intensity bar on the right refers to the latter plot. 
Middle: Doppler map. Contour levels are at 90, 84, and 78 per 
cent of the maximum intensity. Bottom: VI R plot of the origi- 
nal (o) and the reconstructed (x and dashed line) data, and the 
residuals (bottom of the plot). In phases to 1 , error bars are 
given for the original data, phases 1 to 2 show those for the 
reconstructed data. 



Fig. A.IO. Doppler Tomography for AQ Eri, data set Ic. Top: 
Original (left) and difference (original-reconstructed; right) 
spectrum. The intensity bar on the right refers to the latter plot. 
Middle: Doppler map. Contour levels are at 90, 84, and 78 per 
cent of the maximum intensity. Bottom: V/R plot of the origi- 
nal (o) and the reconstructed (x and dashed line) data, and the 
residuals (bottom of the plot). In phases to 1, error bars are 
given for the original data, phases 1 to 2 show those for the 
reconstructed data. 
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Fig.A.ll. Top: Line profiles of selected 10% phase bins for 
VW Hyi. Data sets 1 - 3 are plotted from left to right. Bottom: 
Diagnostic diagram for the combined data set of VW Hyi, al- 
ready corrected for the derived zero phase. The dotted lines 
mark the chosen separation and the corresponding parameters. 
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Fig. A.12. Doppler Tomography for VW Hyi, data set 1. Top: 
Original (left) and difference (original-reconstructed; right) 
spectrum. The intensity bar on the right refers to the latter plot. 
Middle: Doppler map. Contour levels are at 96, 89, and 82 per 
cent of the maximum intensity. Bottom: V/R plot of the origi- 
nal (o) and the reconstructed (x and dashed Une) data, and the 
residuals (bottom of the plot). In phases to 1, error bars are 
given for the original data, phases 1 to 2 show those for the 
reconstructed data. 
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Fig. A.13. Doppler Tomography for VW Hyi, data set 2. Top: 

Original (left) and difference (original-reconstructed; right) 
spectrum. The intensity bar on the right refers to the latter plot. 
Middle: Doppler map. Contour levels are at 97, 90, and 83 per 
cent of the maximum intensity. Bottom: VIR plot of the origi- 
nal (o) and the reconstructed (x and dashed line) data, and the 
residuals (bottom of the plot). In phases to 1, error bars are 
given for the original data, phases 1 to 2 show those for the 
reconstructed data. 



Fig. A.14. Doppler Tomography for VW Hyi, data set 3. Top: 
Original (left) and difference (original-reconstructed; right) 
spectrum. The intensity bar on the right refers to the latter plot. 
Middle: Doppler map. Contour levels are at 90, 80, and 70 per 
cent of the maximum intensity. Bottom: VIR plot of the origi- 
nal (o) and the reconstructed (x and dashed Une) data, and the 
residuals (bottom of the plot). In phases to 1, error bars are 
given for the original data, phases 1 to 2 show those for the 
reconstructed data. 
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Fig. A.15. Top: Line profiles of selected 10% phase bins for RZ 
Leo. Data sets la, lb, 2a, 2b are plotted from left to right. The 
time sequence is symboUsed by the line styles solid - dashed. 
Bottom: Diagnostic diagram of RZ Leo, data set 2, already 
corrected for the derived zero phase. The dotted lines mark the 
chosen separation and the corresponding parameters. 
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Fig. A.16. Doppler Tomography for RZ Leo, data set 1. Top: 
Original (left) and difference (original-reconstructed; right) 
spectrum. The intensity bar on the right refers to the latter plot. 
Middle: Doppler map. Contour levels are at 90, 80, and 70 per 
cent of the maximum intensity. Bottom: V/R plot of the origi- 
nal (o) and the reconstructed (x and dashed Une) data, and the 
residuals (bottom of the plot). In phases to 1, error bars are 
given for the original data, phases 1 to 2 show those for the 
reconstructed data. 
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Fig. A.17. Doppler Tomography for RZ Leo, data set 2. Top: 
Original (left) and difference (original-reconstructed; right) 
spectrum. The intensity bar on the right refers to the latter plot. 
Middle: Doppler map. Contour levels are at 86, 76, and 66 per 
cent of the maximum intensity. Bottom: V/R plot of the origi- 
nal (o) and the reconstructed (x and dashed line) data, and the 
residuals (bottom of the plot). In phases to 1, error bars are 
given for the original data, phases 1 to 2 show those for the 
reconstructed data. 



Fig.A.18. Top: Line profiles of selected 10% phase bins for 
TU Men. Data sets 1 to 4 are plotted from left to right. The 
time sequence is symbolised by the line styles solid - dashed. 
Bottom: Diagnostic diagram for the combined data set of TU 
Men, already corrected for the derived zero phase. The dotted 
lines mark the chosen separation and the corresponding param- 
eters. 
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Fig.A.19. Doppler Tomography for TU Men. Top: Original 
(left) and difference (original-reconstructed; right) spectrum. 
The intensity bar on the right refers to the latter plot. Middle: 
Doppler map. Contour levels are at 94, 84, and 74 per cent of 
the maximum intensity. Bottom: V/R plot of the original (o) 
and the reconstructed (x and dashed Une) data, and the residu- 
als (bottom of the plot). In phases to 1, error bars are given for 
the original data, phases 1 to 2 show those for the reconstructed 
data. 
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Fig.A.20. Top: Line profiles of selected 10% phase bins for 
HS Vir. The time sequence in the left plot is symbolised by 
the sequence solid - dashed. Bottom: Diagnostic for HS Vir, 
already corrected for the derived zero phase. The dotted fines 
mark the chosen separation and the corresponding parameters. 
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Fig.A.21. Doppler Tomography for HS Vir. Top: Original 
(left) and difference (original-reconstructed; right) spectrum. 
The intensity bar on the right refers to the latter plot. Middle: 
Doppler map. Contour levels are at 94, 86, and 78 per cent of 
the maximum intensity. Bottom: V/R plot of the original (o) 
and the reconstructed (x and dashed Une) data, and the residu- 
als (bottom of the plot). In phases to 1, error bars are given for 
the original data, phases 1 to 2 show those for the reconstructed 
data. 



